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In order to identify biomarkers for early diagnosis and/or for therapeutic targets in the delayed health effects 
of ionizing radiation, we analyzed the subgroups of lymphocytes, serum protein levels and gene expression 
profiles in the peripheral blood of three 60 Co y-ray accidentally exposed persons during the three years after ir- 
radiation. Flow cytometry analyses and agarose gel electrophoresis were applied to investigate the subgroups 
of lymphocytes and the composition of serum proteins, respectively. Gene expression profiling was obtained 
using a whole genome gene expression chip assay. Both the percentage of CD4+ T lymphocytes and the ratio 
of Th to Ts were reduced compared with the normal control values. The percentage of albumin decreased 
whereas beta globulin increased. There were 285 up-regulated and 446 down-regulated genes in irradiated 
samples relative to the control samples. The expression of KDR, CEACAM8 and OSM was validated by RT- 
PCR. The majority of the differentially expressed genes encode proteins associated with the immune response, 
inflammation, oncogenesis, cell structure, oxidative stress, neuro-hormone regulation, reproduction, suscepti- 
bility to psychiatric disorders, or transcriptional regulation. We have identified a number of promising novel 
candidates that have potential for serving as biomarkers for delayed damage. Furthermore, the changes in the 
immunological indicator CD4+ T cells, and the ratio of CD4+ T to CD8+ T cells may be biomarkers for the 
prediction of delayed damage by ionizing radiation. The findings of our study are useful for forming a com- 
prehensive understanding of the mechanisms underlying the delayed effects of ionizing radiation. 
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INTRODUCTION 

It is clear that exposure to ionizing radiation can cause 
delayed effects, such as tumorigenesis, cataracts, infertility, 
chronic inflammation, and diffuse fibrosis of the affected 
organ and tissue. Although the human risks of cardiovascu- 
lar diseases and neuro-psychiatric diseases caused by ioniz- 
ing radiation have not been validated by definitive 
epidemiological studies, several reports have described the 
effects of ionizing radiation on the incidence and preva- 
lence of such diseases [1-7]. Understanding cellular and 
molecular pathological changes, and establishing biomar- 
kers for cancer and non-cancer effects of chronic low-dose 
or acute high-dose radiation have become the focus in the 



field of radiation protection in recent years [8-11]. Such 
knowledge will guide the prevention, early diagnosis and 
treatment of ionizing radiation-induced damage, as well as 
the related drug design. It will also provide the theoretical 
basis for adopting efficient measures in the practice of radi- 
ation protection (e.g. health surveillance of radiation 
workers), medical treatment and therapeutic strategies in ra- 
diation emergency medicine, and for providing health ser- 
vices for long-term cancer survivors after radiotherapy. 

According to the National Cancer Prevention office, there 
were 2 million new cancer occurrences in the year 2000 in 
China. Among them, 400 000 cases received radiotherapy. 
After secondary malignancies, cardiovascular disease is the 
leading cause of late morbidity and death among cancer 
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survivors in the United States [12]. Similarly, we can specu- 
late that the number of cancer survivors is increasing in 
China, even though there is no up-to-date statistical data 
available, and that non-cancer diseases will occupy an im- 
portant role in late morbidity and death in cancer patients 
after radiotherapy and chemotherapy. From 1994-1998 in 
China there were 139 cases of radiation accidents and 238 
persons were overexposed to radiation to some degree. There 
were 230 000 radiation workers in China up to the year 
2009. Among them, 1718 were identified as sustaining 'radi- 
ation damage', and 464 were diagnosed as having 'suspected 
radiation disease'. Therefore, there is a great need to investi- 
gate delayed radiation effects in the above-mentioned three 
categories of persons, particularly since the use of the emer- 
ging tools of bioinformatics technology has shown that the 
effects of large doses are distinguishable from those of low 
doses. The early identification of risk groups and the diagno- 
sis of oncological and degenerative diseases will be highly 
beneficial. Therefore, we performed analyses using blood 
samples from three persons who were accidentally exposed 
to high-dose radiation, simulating the biological condition of 
cancer survivors who have been treated with high-dose radi- 
ation, and those who have been accidentally exposed to radi- 
ation in some other way. 

The three accidentally exposed persons followed up in 
the present study were exposed in a radiation-sterilization 
facility accident that occurred on 12 May 2008. Five male 
workers from 38^5 years of age were exposed to Co-60 
source (activity, 1.6415 x 10 4 Ci) y- ra y s f° r nearly 20 
minutes. The Co-60 source was being used for sterilization 
of Chinese herbs. The initial cause of this accident was a 
loss of function of the facility safety interlock. The five 
workers entered the radiation room to move out the irra- 
diated herbs and move in new herbs for irradiation, while 
failing to check the radiation meter count. Additionally, the 
audio alarm sound of the portable radiation meter was 
masked by the noise of the ventilation system in the radi- 
ation room. Only when one person noticed that the steel 
line controlling the up-and-down movement of the radiation 
source was very loose did they become aware that the 
source was still in exposure status. They evacuated the site 
after they dropped the source down to the safe position. 
The five persons were sent to the Affiliated Hospital of 
China Institute for Radiation Protection 40 min later. 
Physical doses were estimated using the thermolumines- 
cence dosimetry system in the Dosimetry Laboratory of the 
institute, according to the simulated exposure scenario. The 
physical doses were 12.0, 11.2, 1.82, 1.82 and 1.0 Gy, re- 
spectively. The five persons were sent to the Affiliated 
Hospital of Beijing Institute of Radiation Medicine for 
further treatment. They were diagnosed with acute radiation 
syndrome, one with gastrointestinal syndrome, and the 
other four with hematopoietic syndrome with either a high 
or medium level of severity. Four of the five persons 



received bone marrow transplantation from their brothers or 
son. Two persons with physical doses of 12.0 and 1.82 Gy 
died on Day 60 and 17 months after the accident, 
respectively. 

The three remaining persons (here named A, B and C) 
were followed annually with clinical examinations, labora- 
tory tests and biodosimetric estimations. Among them, 
A and B received bone marrow transplants on Day 14 and 
21 after the accident, respectively. Patient A presented with 
the most serious injuries, manifesting acute double pneu- 
monia, upper respiratory infection and middle ear infection 
during the treatment. His lymphocyte percentage went back 
to the normal level more than two years after the accident. 
Patient B experienced acute pneumonia and upper respira- 
tory infection. His lymphocyte percentage went back to the 
normal level within six months. Patient C was treated with 
G-CSF 48 h after the accident and manifested upper re- 
spiratory infection within two months. The three patients 
were monitored for delayed health effects. We noticed 
some diseases appeared gradually in the exposed persons, 
including opacity of the lens, skin disorders, sperm abnor- 
malities, and abnormalities in lymphocyte subgroups, as 
well as osteopenia. The mechanisms underlying these 
delayed effects of exposure to a high dose of irradiation are 
not well known. Previous studies demonstrated that oxida- 
tive stress is associated with the acute skin reactions caused 
by high doses of irradiation [13-15], and changes in 
immunological indicators were found in the 'liquidators' at 
the Chernobyl accident and in the atomic bomb survivors 
[1-3]. Therefore, the aim of the present study was to under- 
take follow-up studies of the immunological status and 
gene expression profiles of the three accidentally exposed 
persons in order to establish candidate biomarkers for the 
early identification of persons at risk or for the diagnosis of 
oncological diseases and other delayed harmful effects 
caused by ionizing radiation. 

MATERIALS AND METHODS 

Reagents 

Cy3 NHS ester was from GE, Healthcare (Connecticut, 
USA) and aaUTP was from Ambion (California, USA). The 
fluorochrome-conjugated antibodies CD4-FITC/CD8-PE/ 
CD3-PC5, CD3-FITC/CD(16 + 56)-PE and CD19-PC5, and 
the isotypic antibodies for control IgGl-FITC/IgGl-PE/ 
IgGl-PC5 were from Immunotech, Beckman, France. 
Ethanol and P-diathiothretol were from Sigma. RNase-free 
water, and Nuclease-free water were products of HyClone, 
(Logan, UT, USA). Water was double distilled and deionized 
in an ultra-pure water apparatus (Millipore). 

Blood sampling and preparation 

Fresh blood samples were collected from the three persons ac- 
cidentally exposed to 60 Co-y-rays (three years post-irradiation) 
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and three age-matched, non-irradiated healthy donors. 
Leukocytes were isolated and total RNA was extracted. 
Four RNA samples, one from each of the three accidentally 
exposed persons, along with a RNA mixture extracted from 
a combination of three healthy donors' leukocytes (ratio 
1:1:1, each 3.3 x 10 6 leukocytes), were subjected to Agilent 
4 x 44 k whole genome gene expression microarray assay. 
Fresh blood samples collected from the three exposed 
persons were used for the leukocyte count, measurement of 
the subgroups of lymphocytes by flow cytometry, and 
separation and quantification of serum proteins by agarose 
gel electrophoresis. The measurements were conducted in 
The Second Affiliated Hospital of Shan Xi Medical 
College. The experimental systems with the control refer- 
ence levels that we used in the present study are routinely 
used for clinical diagnosis in the hospital. All the assays 
were performed following a single blood draw during 
yearly follow up study. The overview of the experimental 
design is shown in Fig. 1. 



Whole genome gene expression microarray assay 

Total RNA was extracted using the QIAamp RNA Blood 
Minikit (QIAGEN GmbH, Hilden, Germany) following the 
manufacturer's instructions. An RNase-Free DNase Set 
(QIAGEN) was used to digest DNA efficiently during 
RNA purification. RNA quality was assessed by spectro- 
photometry (NanoDrop, ND1000, Massachusetts, USA) 
and gel electrophoresis. To be considered high quality, the 
RNA had to have a 260/280 ratio > 1.7 and an 18S/28S 
rRNA ratio ~2. 

Labeling was done using a Low RNA Input Linear 
Amplification Kit (Agilent, California, USA) according to 
the manufacturer's instructions, followed by RNAeasy 
Column Purification (QIAGEN, Hilden, Germany). Briefly, 
cDNA was generated using a T7 promoter element coupled 
oligo dT primer. From the cDNA, labeled cRNA was gen- 
erated via an in vitro transcription reaction using T7 RNA 
polymerase and aa-UTR aa-UTP incorporated into cRNA 
was dyed with Cy3 NHS ester. Dye incorporation and 
cRNA yield were checked with the NanoDrop ND-1000 
Spectrophotometer. 



Labeled cRNA (825 ng) from the respective sample was 
fragmented and hybridized onto Agilent whole human 
genome microarray s (4x44K G4112F, slide number, 
66634, California, USA ) using a Gene Expression 
Hybridization Kit (5188-5242, Agilent Technologies). 
Hybridization was performed for 17 h, rotating at a speed 
of 10 rpm at 65 °C in an Agilent hybridization oven 
(G2545A). After hybridization, microarray s were washed 
for 1 min at room temperature with GE Wash Buffer 1, and 
1 min with 37°C GE Wash buffer 2 (5188-5327, Agilent 
Technologies). 

The slides were scanned with Agilent microarray scanner 
(G2505C) using one color setting for 4 x 44 K array slides 
(scan resolution 5 um, dye channel set to green, and green 
PMT set to 100%). 

Images were processed with Agilent Feature Extraction 
Version 10.7.1.1 to obtain raw data files. The raw data files 
that contained 'Processed Signal' were processed using the 
Limma package developed within the Bioconductor project 
in the R statistical programming environment. After log 
transformation, the data were normalized using the 'quan- 
tile' method. The gene expression data were filtered accord- 
ing to the values of AFE quantification flags of glsFound, 
IsWellAboveBG, glsSaturated and the computing results of 
WellAboveNeg. After combinations of the replicated 
probes, the processed data was used to analyze the differen- 
tially expressed genes using a threshold of absolute log2 
fold change value cut off at 1 . 



RT-PCR analysis 

The quantitative RT-PCR was carried out with a 7500 
Real-Time PCR System (Applied Biosystems, California, 
USA), with primers for KDR, OSM, CEACAM8 and 
ACTB, as shown in Table 1. The primers were designed 
using Primer Express 2.0. The synthesis of cDNA and 
fluorescence PCR amplification were performed using an 
one-step QuantiFast SYBR Green RT-PCR Kit (QIAGEN). 
Relative quantification of mRNA expression was performed 
by the delta delta Ct (CtO) method using the level of 
P-actin mRNA for normalization of mRNA abundances 
across the samples. 2~ Ct( - Ct(y> was used for calculation of the 
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Table 1. 


Primers for genes of KDR, OSM, CEACAM8 and ACTB 




Gene 


Sequence 


Length of product (bp) 


KDR 


forward: 5 ' CTC AAACCAGACA AGCGGCTA3 ' 
backward: 5 ' GCTACCGGTTTGCACTCC AAT3 ' 


121 


OSM 


forward: 5'CCAGACTTCCTCCTTTCCGTG3' 
backward: 5 ' AC ACCCTGCCGCTGTTACAG3 ' 


137 


CEACAM8 


forward: 5 ' TC AGTGACCCAGTCACCCTGA3 ' 
backward: 5 ' TGG ATTAG AGGCCGC ATGG3 ' 


122 


ACTB 


forward: 5 ' A AGTACTCCGTGTGG ATCGGC3 ' 
backward: 5 ' GCCTAGAAGCATTTGCGGTG3 ' 


125 



relative expression amount: 

CtO Ctt ar g e t gene Cti n t e mal standard gene ; 

Ct Ct tes t Ct conu - 0 i 



Measurement of serum proteins by agarose gel 
electrophoresis 

The serum proteins were separated and quantified by 
agarose gel electrophoresis using the REP SP Applicator 
Kit (Helena, Biosciences Europe). After programming the 
parameters into the automatic instrument Spife Combo 
(Helena Laboratories, USA), 40 ul of the control/sample 
serum was pipetted into the appropriate sample well and 
placed in the correct position on the sample tray. Following 
electrophoresis, the gel was stained with acid blue stain for 
10 min, destained with 5% glacial acetic acid and dried at 
65°C in a drying oven. The REP Applicator kit separates 
proteins into five main classes (albumin, alpha 1 -globulin, 
alpha2-globulin, beta-globulin and gamma globulin) 
according to charge in an agarose gel. The stained gel was 
subjected to visual and quantitative interpretation. 

Measurement of subgroups of lymphocytes by flow 
cytometry 

Three reaction tubes were designated for every blood 
sample, and contained OptiClone monoclonal antibodies 
CD4-FITC/CD8-PE/CD3 -PC5 , CD3-FITC/CD(16 + 56)-PE 
and CD19-PC5, and corresponding isotypic controls, re- 
spectively. The 12x75 mm tubes contained 20 ul of 
OptiClone monoclonal antibodies, and 100 ul of anticoagu- 
lated whole blood was added and briefly vortexed to mix. 
The reactions were incubated for 15 min at room tempera- 
ture (18-25°C), in the dark. The red blood cells were lysed 
using Optilyse lysing solution, and incubated for 10 min at 
room temperature, protected from light. Following centrifu- 
gation for 5 min at 150 xg at room temperature, the cell 
pellets were washed twice with 3 ml of PBS. The resulting 



cell pellets were dissolved in 0.5 ml PBS. Acquisition and 
analysis were performed with a Cytometric FC 500 
(Beckman Coulter, California, USA) flow cytometer 
equipped with CXP Software. 

RESULTS 

Changes in the subgroups of lymphocytes 

The subgroups of lymphocytes in the peripheral blood of 
the three exposed persons changed during the three years 
following the accident. The leukocyte count was elevated 
in all three individuals within hours after the exposure but 
returned to normal after one year. The percentage of lym- 
phocytes was depressed in all three individuals within 
hours after the exposure and remained depressed in one in- 
dividual for two years (Table 2). Changes in lymphocyte 
subgroups are shown in Table 3. The percentages of helper 
T (Th, CD3 + CD4+) cells were reduced in all three 
persons. The ratios of Th to Ts (CD3 + CD8 + , Suppressor 
T) were also reduced compared with the reference level. 
The total T cells decreased or were near the lower limit of 
the reference level. 

Changes in serum protein levels 

The percentage of albumin was lower, whereas Beta globu- 
lin was higher compared with the reference values 
(Table 4). 

Differentially expressed genes 

There were 731 differentially expressed genes in the exposed 
persons. Among these genes, 285 were up-regulated and 446 
were down-regulated. The raw data had been accepted by 
the GEO database of NCBI (http://www.ncbi.nlm.nih.gov/ 
geo/ Access number GSE34131). Tables 5-7 show some of 
the up-regulated, down-regulated and immunological func- 
tion-associated genes, along with protein-associated informa- 
tion including subcellular localization, biological processes, 
molecular function and/or association with diseases (see 
Supplementary Materials). The absolute log2 ratio of the 
gene lists was cut off at 2.0. 
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Table 2. Leukocytes and lymphocyte percentage in the peripheral blood of accidentally exposed persons 



Cell type Time after radiation 


A 

(11.2 Gy) 


B 

(1.82 Gy) 


C 

(1.0 Gy) 


Reference value 


Leukocyte (xl0 9 /l) 


2h 


9.89 


7.97 


10.4 | 


4-10 10 9 /1 




4 h 


12.75 t 


12.09 t 


15.86 T 






1 year 


5.73 


4.61 


5.01 






2 year 


7.19 


4.4 


5.52 






3 year 


5.46 


4.67 


5.57 




Lymphocyte (%) 


2h 


8.921 


14.621 


12.821 


20-40% 




4 h 


3.52 J. 


4.821 


4.62 | 






1 year 


18.02 1 


25.8 


25.7 






2 year 


16.02 1 


27 


24.5 




Table 3. Lymphocyte sub-groups distribution (%) in 


the accidentally exposed persons 






Cell classification 


Time after 
irradiation (year) 


A (11.2 Gy) 


B (1.82 Gy) 


C (1.0 Gy) 


Reference value (%) 


Th (CD3 + CD4 + ) 


1 


18.7| 


25.91 


26.7| 


30.09-40.41 




2 


19.0| 


27.11 


28.0| 






3 


19.51 


28.2 X 


27.2| 




Ts (CD3 + CD8 + ) 


1 


41.7 t 


27.8 


29.6 


20.74-29.42 




2 


34.6 T 


25.5 


30.9 t 






3 


35.9 T 


24.3 


29.3 




Th/Ts 


1 


0.45 | 


0.931 


0.90| 


0.98-1.94 




2 


0.55| 


1.06 


0.91 | 






3 


0.541 


1.16 


0.93 | 




Total T cell (CD3+) 


1 


62.1 


57.61 


59.9| 


60.0-80.0 




2 


56.1 I 


56.4 J. 


61.8 






3 


56.61 


56.0 J. 


60.1 




NK cell {CD3-CD(16 + 56) + ) 


1 


5.501 


24.5 


27.8 t 


8.1-25.6 




2 


5.0 { 


28.3 T 


22.2 






3 


5.3 I 


23.6 


17.2 




CD3 + CD(16 + 56) + 


1 


0.30 


0.10 


0.30 


<5 




2 


8.6 T 


5.4 T 


3.4 






3 


4.8 


3.0 


1.4 




CD19 + 


1 


24.6 T 


9.40 


13.9 


9-15 




2 


38.1 T 


14.8 


15.1 






3 


34.5 T 


16.6 T 


11.5 





RT-PCR analysis 

The differentially expressed genes CEACAM8, OSM and 
KDR were validated by real time RT-PCR analysis. The ex- 
pression analysis showed an up-regulation of CEACAM8 



and OSM, and a down-regulation of KDR, relative to 
control levels (Supplementary Table 8). Thus the expres- 
sions of the three genes were consistent with that of whole 
genome gene expression CHIP assay. 
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Table 4. Serum proteins distribution (' 


%) in the three accidentally exposed persons 






Component 


Time after 
irradiation (year) 


A (11.2 Gy) 


B (1.82 Gy) 


C (1.0 Gy) 


Reference value (%) 


Albumin 


1 


58.65 


46.891 


f\l 1A 
OZ. / 


^n f\f\ nn 

JZ. J\J— OD.UU 




2 


42.751 


50.891 


42.261 






3 


49.491 


45.591 


49.251 




Alpha 1 


1 


3.45 


6.43 


1 Q7 1 

Z.V 1 i 


JJU-/.UU 




2 


5.06 


2.551 


4.36 






3 


5.99 


4.41 


2.551 




Alpha2 


1 


7.59 


7.23 




rS ^0 1 1 70 

U . J U 1 1 . / U 




2 


15.711 


8.78 


10.79 






3 


10.42 


6.83 


9.19 




Beta 


1 


16.041 


12.58 


19 81 






2 


18.261 


14.831 


16.441 






3 


17.371 


16.271 


17.061 




Gamma 


1 


14.27 


26.871 


1 f\ 11 

i o. / z 


11 1 n in ah 




2 


18.21 


22.96 


26.16| 






3 


16.73 


26.901 


21.49T 




Ratio 


1 


1.42 


0.88 


1.68 






2 


0.75 


1.04 


0.73 






3 


0.98 


0.84 


0.97 





DISCUSSION 

There are few validated biomarkers for the prediction or 
early diagnosis of delayed health effects caused by expos- 
ure to ionizing radiation. The knowledge of characteristic 
alterations occurring at the level of the tissue, cell and mol- 
ecule after radiation will help identify the pathological pro- 
cesses involved in delayed harmful effects. Radiation 
biomarkers are indicative components that play important 
roles in the signaling pathways underlying the pathological 
processes. There are few reports in the literature on radi- 
ation biomarkers in accidentally exposed persons. The 
goals for identifying radiation biomarkers are: (i) to charac- 
terize the mechanisms underlying delayed radiation 
damage; (ii) early diagnosis or prediction of delayed 
effects; (iii) assessment of the risk of secondary health 
effects in cancer patients treated with radiotherapy. 
Generally, radiation biomarkers are categorized by applica- 
tion into indicators for diagnosis, dose estimation, sensitiv- 
ity and risk assessment. 



In the present study we followed biomarkers for three 
years in three accidentally exposed persons. Our results 
showed that the percentage of CD4+ T lymphocytes 
remained lower and the ratio of CD4+ T lymphocytes to 
CD8+ T lymphocytes decreased compared with the refer- 
ence levels. The level of albumin (known to be involved in 
the non-specific immune response) decreased, whereas the 
level of beta globulin increased. These changes were per- 
sistent for the three years after the accident. Further, gene 
expression profiles in the leukocytes indicated that the dif- 
ferentially expressed genes were involved in the immune re- 
sponse, inflammation, oncogenesis, cell structure, oxidative 
stress, neuro-hormone regulation, reproduction, susceptibil- 
ity to psychiatric diseases, and regulation of gene expres- 
sion, which were summarized by reviewing the 
Geneontology, PubmedOMIM and Pubmedgene databases 
and related original research papers. Down-regulated genes 
predominated among the differentially expressed genes. 

Epidemiological analysis of the immune status in the 
'liquidators' of the Chernobyl accident showed that the 
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Table 5. Partial list of uprej 


julated genes 






/\l-l,c»MUll IIUIHUCI 


vrCIlG 


11 U LI 11 1 




MM 01870? 

IN 1V1 UIO / Ui. 


ADARR? 

r\ YJ t\ I\ I J Z. 


nHi^TinciTii 3 np'itnin'jcp R9 
dU-OllUolllt U.CdiilllldftC DZ 


u 


MM 001 09^ 
IN 1V1 UU 1 yL J 


TYPFA A 
UCiTn'r 


defensin, alpha 4 


J 


NM 006671 


CT pi A 7 


c r\l 1 1 tt* c i t~t~i #^»T" fo in il\7 1 i cYln t Qm ctt^* tret n c T"\^t"f i mpmnpr / 
SUILILC UdlilCl Idllllly 1 ^glULdllldLC LidllspUl LCl^, lilClllDCl / 


J.J 


MM 004 1 7 1 

IN 1V1 UUH- 1/1 


<iT P1 A 9 


solute carrier family 1 , member 2 


0 7 


mm oooss^ 


VjiJ 111 


glutathione S-transferase theta 1 


4 7 


MM 007,19 A 
IN 1V1 UUZH-ZH- 


MMP8 
JVllvlr o 


matrix metallopeptidase 8 


j.y 


mm nrw 1 A'x 


Tip 
ill 


haptoglobin 


A 1 
4-.J 


MM 09000^ 
In lvl_UZUy y D 


I TOO 

I I IK 


haptoglobin-related protein 




MM OO'SQRO 

IN 1V1 UUJ70U 


SI OOP 


>I 1 00 r*nlr*inm nitinincr nrntpin P 
OlUU Ud.lClU.lll L/llllllllii U1ULC111 IT 


3 7 
j. / 


MM 00 1 Q79 
IN 1V1 UU 1 y 1 L 


"PT AMP 


elastase 


Z.J 


mm 0004Q4 


TOT 17A1 


rrtllnopn tvnp ^VTT nlnhn 1 
V/Uiiugcii, l y ut yv v ±±5 ciii 'iiii j. 


3.0 


NM_004496 


F0XA1 


forkhead box Al 


3.9 


NM_015714 


G0S2 


G0/G1 switch 2 


3.7 


NM_004430 


EGR3 


early growth response 3 


3.0 


NM_000963 


PTGS2 


inducible cyclooxygenase 


2.5 


NM_001549 


IFIT3 


interferon-induced protein with tetra-tricopeptide repeats 3 


2.0 


Table 6. Partial list of down-regulated genes 


Accession number 


Gene 


Protein 


Log2Ratio 


NM_198129 


LAMA3 


laminin, alpha 3 


-4.3 


NM_001464 


ADAM2 


ADAM metallopeptidase domain 2 


-4.7 


NM_001079858 


GPR64 


G protein-coupled receptor 64 


-4.6 


NM_194250 


ZNF804A 


zinc finger protein 804A 


-4.1 


NM_ 145740 


GSTA1 


glutathione S-transferase alpha 1 


-4.0 


NM_0 16931 


NOX4 


NADPH oxidase 4 


-4.3 


NM_001 137669 


RGSL1 


regulator of G-protein signaling like 1 


-4.1 


NM_005438 


FOSL1 


FOS-like antigen 1 


-3.4 


NM_020868 


DPP 10 


dipeptidyl peptidase 10 


-3.9 


NM_000817 


GAD1 


glutamate decarboxylase 1 


-3.4 


NM_153699 


GSTA5 


glutathione S-transferase alpha 5 


-3.9 


NM_004883 


NRG2 


neuregulin2 


-3.4 


NM_003551 


NME5 


non-metastatic cells 5 


-3.2 


NM_001836 


CMA1 


chymase 1 


-3.4 


NM_00 1040 105 


MUC17 


mucin 17 


-3.4 


NM_080801 


C0L13A1 


collagen, type XIII, alpha 1 


-3.0 


NM_002538 


0CLN 


occludin 


-3.4 



percentage of CD4+ T lymphocytes decreased, and the 
ratio of CD4+ to CD8+ T lymphocytes decreased; the abso- 
lute number of CD8+ lymphocytes increased over one to 



three years before a cancer-confirming diagnosis [8]. 
Investigation of the immune status and inflammatory re- 
sponse in the atomic bomb survivors showed that the 
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Table 7. Partial list of immunological function-associated genes 



Accession number 


Gene 


Protein 


Log2Ratio 


NM 000588 


IL3 


interleukin 3 


-3 7 


RC040272 


IL16 


lntprlpukin 1 f\ 


-2 9 1 


NM_002253 


KDR 


kinase insert domain receptor 


-3.5 | 


NM_002175 


IFNA21 


interferon, alpha 21 


-2.5 | 


NM_030788 


TM7SF4 


transmembrane 7 superfamily member 4 


+2.5 T 


NM_020530 


OSM 


oncostatin M 


+2.6 T 


NM_001816 


CEACAM8 


carcinoembryonic antigen-related cell adhesion molecule 8 


+2.5 T 


NM_002343 


LTF 


lactotransferrin 


+2.5 T 


NM_000547 


TPO 


thyroid peroxidase 


+2.8 T 


NM_080657 


RSAD2 


radical S-adenosyl methionine domain containing 2 


+2.6 T 



immunity of the T cells was compromised, including an 
impaired mitogen-stimulated proliferation, a decrease in the 
quantity of naive T cells, an increase in the function-weak 
memory CD4+ T cells, and an increase in CD25 + / 
CD 127- Reg T (Regulatory T) cells in a dose-dependent 
way [9]. The roles of CD4+ and CD8+ T cells, 
B-lymphocytes and natural killer cells in the prediction of 
radiation-induced delayed toxicity in cervical cancer have 
been reported [10]. That research data, combined with our 
results showing alterations in subgroups of lymphocytes, 
serum proteins and gene expression in the accidentally 
exposed persons, suggests that persistent immunological 
change is one of the main processes after a high dose of ra- 
diation, and it may also be an important focus for revealing 
the mechanisms underlying the delayed effects of ionizing 
radiation. Moreover, the decrease in CD4+ T cells and the 
decrease in the ratio of CD4+ to CD8+ T cells may be bio- 
markers for prediction of delayed damage from ionizing 
radiation. 

Our results from the agarose gel electrophoresis analysis 
showed that beta globulin increased. Lactotransferrin (LTF) 
(Table 7) and haptoglobin (HP) (Table 5), two components 
of the beta globulin band area, are 2.5 and 4.0 times higher 
in mRNA levels relative to the control levels. This suggests 
that the increase in LTF and HP production may contribute 
to the increase in protein production in the beta band area, 
indicative of activation of the non-specific immune re- 
sponse and inflammatory response. 

Except for genes EGR3, IL3, IL16, KDR, GSTT1, AD ARB 2, 
and PTGS2 that have been reported as radiation-associated 
[16-21], most of the genes shown in Tables 5-7 are shown for 
the first time to be involved in radiation-induced responses. 

113, IL16 and KDR are involved in the regulation of the 
proliferation, growth and differentiation of hematopoietic 
cells, T cells and endothelial cells. The down-regulation of 
those genes, as shown in Table 7, suggests a cell-growth 
inhibitory effect has been induced by radiation. Similarly, 



the down-regulation of NRG2 and FOSL1 (shown in 
Table 6) reflects the inhibitory effect on cell growth. NRG2 
induces growth and differentiation of epithelial, neuronal 
and glial cells [22]. FOSL is heterodimerized with proteins 
of the JUN family forming transcription factor complex 
AP-1, and is implicated as a regulator of cell proliferation, 
differentiation and transformation [23]. The up-regulation of 
EGR3 (shown in Table 5) may present negative or positive 
regulation of cell growth and differentiation [17]. 

CEACAM8 is known as the activation marker of human 
granulocytes and regulates adhesion and activation of 
human eosinophils [24]. TM7SF4, preferentially expressed 
by dendritic cells, presents as an antigen to naive T cells 
[25]. OSM and OSMR function to regulate cytokine produc- 
tion (including IL-6, G-CSF and GM-CSF) by endothelial 
cells, and are also known as inflammatory markers [26]. 
Thus, the up-regulation of those genes in our results sug- 
gests activation of granulocytes, T cells and immune 
responses in the accidentally exposed persons. 

The up-regulation of HP, HPR, PTGS2, IFIT3, OSM, 
OSMR and MMP8, which is associated with inflammation, 
reflects the occurrence of the inflammatory response in the 
accidentally exposed persons. 

The alterations to the glutathione metabolic pathway are 
intriguing. The function of the glutathione-s-transferase 
family is to catalyze the conjugation of reduced glu- 
tathiones and a variety of electrophiles, such as mutagens 
and carcinogens, thereby playing roles in detoxification. 
The changes in expression levels of these family members 
in the accidentally exposed persons might be attributed to 
effects resulting in perturbations in oxidative metabolism. 
Such effects may lead to excess generation of reactive 
oxygen species, including free radicals. The latter effect 
may be associated with altered expression or activity of 
antioxidants. Among the family members, GSTT1 was 
up-regulated whereas GSTA1 and GSTA5 were down- 
regulated. One study showed that deletion of GSTT1 
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increased the risk of senile cataracts in Asian people [27]. 
The up-regulation of GSTT1 in the accidentally exposed 
persons might be accounted for as a functional compensatory 
increase to protect the body from forming cataracts, and this 
is supported by the occurrences of opacity of the lens in the 
three accidentally exposed persons at different times after the 
accident. Glutamate transporters of SLC1A7 and SLC1A2 
were up-regulated, and they might serve to clear the excita- 
tory neurotransmitter glutamate from the extracellular space 
of synapses in the central nervous system [28]. GAD1 was 
down-regulated in our results, and according to the literature 
[29, 30], down-regulation of this gene has been seen in 
Parkinson's disease and schizophrenia. 

The up-regulation of ADARB2, a RNA editing enzyme, 
and FOXA1, a transcriptional activator for liver- specific tran- 
scripts such as albumin and transthyretin, suggests enhance- 
ment of the expression of genes that regulate a multitude of 
responses, which when perturbed leads to tumorigenesis. 

CEACAM8 was up-regulated; its high expression has been 
found in acute lymphoblastic leukemias [31]. MUC17 was 
down-regulated; its low expression has been found in inflam- 
mation and neoplastic disease in the digestive tract [32]. 
RGSL1 was down-regulated, and its function is to increase 
the GTPase activity of the G a subunit to attenuate signaling 
from the G-protein-coupled receptor, located in lq25.3, the 
prostate cancer region (HPC1) associated with breast cancer 
and prostate cancer [33]. SI OOP was up-regulated; its high 
expression has been reported to be associated with prostate 
cancer progression [34]. A high level of expression has also 
been found in pancreatic ductal adenocarcinoma and SWOP 
may be a marker for early diagnosis [35]. OCLN was down- 
regulated; it encodes an integral membrane protein required 
for the cytokine-induced regulation of the tight junction 
paracellular permeability barrier, and one report has shown 
that loss of OCLN leads to the progression of human breast 
cancer [36]. 

GPR64, NME5 and ADAM2 were down-regulated. 
GPR64 is an epididymis-specific transmembrane protein, 
NME5 exhibits a nucleoside diphosphate kinase activity 
[37] and is involved in spermiogenesis and flagellar move- 
ment [38], and ADAM2 is a subunit of an integral sperm 
membrane glycoprotein called fertilin that plays an import- 
ant role in sperm-egg interactions. The down-regulation of 
those three genes was shown to coincide with the occur- 
rences of poor sperm motility, low sperm count, or semen 
with no sperm, as shown by the clinical examinations of 
the accidentally exposed persons. 

DEFA4 was up-regulated in our results. It is found in the 
granules of neutrophils and has corticostatic activity, inhi- 
biting corticotrophin-stimulated corticosterone production. 
It functions via the renin-angiotensin system pathway. 
CMA1 was down-regulated in our results. It is expressed in 
mast cells, and is largely responsible for converting angio- 
tensin to the vasoactive peptide angiotensin. Angiotensin II 



has been implicated in blood pressure control, and a per- 
turbation of its yard is associated with the pathogenesis of 
hypertension, cardiac hypertrophy, and heart failure. The 
changes in expression level of these two genes suggest a 
response of the renin-antiotensin system to ionizing radi- 
ation in the accidentally exposed persons. SLCA2 was 
up-regulated and its encoded protein functions to clear the 
excitatory neurotransmitter glutamate from the extracellular 
space at synapses in the central nervous system. Glutamate 
clearance is necessary for proper synaptic activation and to 
prevent neuronal damage from excessive activation of glu- 
tamate receptors. Decreased expression of this protein or 
loss of its function is associated with amyotrophic lateral 
sclerosis, Parkinson's disease, Alzheimer's disease, autism 
spectrum disorders and psychiatric disorders [28]. Taken to- 
gether, those data suggest the disturbance of the neuro- 
hormone regulation system and a susceptibility to neuro- 
psychiatric disorders in the accidentally exposed persons. 

LAM A3 was down-regulated in our results. Its encoded 
protein is a component of basement membrane, involved in 
cell adhesion, signal transduction and differentiation of kerati- 
nocytes. COL17A1 was up-regulated. It encodes a transmem- 
brane protein, a structural component of hemidesmosomes, 
multiprotein complexes at the dermal-epidermal basement 
membrane zone that mediate adhesion of keratinocytes to the 
underlying membrane. Mutations in LAMA3 and COL17A1 
have been associated with junctional epidermolysis bullosa 
[39, 40]. Skin disorders manifesting pigmentation in the face, 
desquamation and itching in the skin of the arms and lower 
legs, thickening of the fingernails, or psoriasis that appeared 
consecutively after the exposure accident might be related the 
alterations of those two genes. 

We have discussed the immune response, the inflamma- 
tory response, oxidative stress, tumorigenesis, spermiogen- 
esis abnormality, skin disorders and disturbance of 
neuro-hormone regulation axes in the accidentally exposed 
persons in the above paragraphs. Additionally, our results 
suggest that CD4+ T cells, and the ratio of CD4+ to CD8+ 
T cells may be biomarkers for radiation-caused delayed 
damage. Collectively, we propose the following candidate 
genes as biomarkers for delayed damage: (i) inflammation 
and immunity associated genes: IL3, IL6, HP, HPR, 
PTGS2, IFIT3, ELANE, OSM, OSMR and MMP8; (ii) 
oncogenesis associated genes: CEACAM8, SIOOP, MUC17, 
RGSL1 and OCLN; (iii) the glutathione-s-transferase family 
and the metabolism of glutamate associated genes: GSST1, 
GSTA1, GSTA5, SLC1A2, SLC1A7 and GAD1; (iv) sperm- 
atogenesis associated genes: GPR64, NME5 and ADAM2; 

(v) skin disease associated genes: LAMA3 and COL17A1; 

(vi) neuro-hormone regulation and susceptibility of psychi- 
atric disorders associated genes: DEFA4, CMA1, GAD1 
and SLCA2. The alterations of expression of the above 
genes eventually lead to the onset of delayed injuries. 
Therefore, those genes are candidates for the prediction of 
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Fig. 2. Differentially expressed genes, related pathological processes and delayed effects of ionizing radiation. The diseases (late effects) 
in the shaded rectangles indicate the diseases already occurring in the exposed individuals, while the diseases noted in black ellipses are 
as yet non-occurring. Solid arrows show the results of gene expressions, and dotted arrows connect the pathological processes indicated 
by the results of gene expressions to the delayed effects (diseases). 



delayed injuries. The relationship between differentially 
expressed genes, related pathological processes, and 
delayed effects (diseases) is outlined in Fig. 2. 

One issue we should note is that two of the three acci- 
dentally exposed persons were given bone marrow trans- 
plantation, the other one was administered with G-CSF. 
The influence of those treatments on the delayed effects is 
unknown. It is clear, however, that those appropriate treat- 
ments improve survival in LD50/30 exposed persons. 
Moreover, it should be recognized that the gene expression 
chip assay was not repeated for confirmatory purposes and 
thus generation of an indicator of statistical significance has 
not been provided. 

The candidate genes noted above are members of signal- 
ing pathways in processes responding to a high dose of 
gamma rays. Their biological functions and association 
with diseases will become clearer with increasing time. 
This will ultimately allow characterization of biomarkers 
that will permit prediction or early diagnosis of delayed 
harmful effects of radiation to be established more 
accurately. 

SUPPLEMENTARY DATA 

Supplementary data is available at the Journal of Radiation 
Research online. 
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